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SCALING LIMITS FOR CONTINUOUS OPINION 
DYNAMICS SYSTEMS 

By Giacomo Como and Fabio Fagnani 

Massachusetts Institute of Technology and Politecnico di Torino 

Scaling limits are analyzed for stochastic continuous opinion dy- 
namics systems, also known as gossip models. In such models, agents 
update their vector-valued opinion to a convex combination (possi- 
bly agent- and opinion-dependent) of their current value and that 
of another observed agent. It is shown that, in the limit of large 
agent population size, the empirical opinion density concentrates, at 
an exponential probability rate, around the solution of a probability- 
measure- valued ordinary differential equation describing the system's 
mean-field dynamics. Properties of the associated initial value prob- 
lem are studied. The asymptotic behavior of the solution is analyzed 
for bounded-confidence opinion dynamics, and in the presence of an 
heterogeneous influential environment. 

1. Introduction. In this paper, we undertake a rigorous mathematical 
analysis of a family of stochastic dynamical systems proposed as opinion dy- 
namics models in the recent literature: see, for example, [11], Section III, [24], 
and references therein. Here, we shall focus on the so-called "gossip" mod- 
els, where the information propagation, as the name suggests, takes place 
through pairwise interactions. These models have been proposed in other 
scientific areas, for instance, as aggregation and estimation algorithms in 
sensor and robotic networks (see, e.g., [9, 28]), or as models for aggregation 
and clustering in biological systems (see, e.g., [16]). 

One of the simplest gossip model can be described as follows. Each agent a 
of a population A of finite size n := |^| possesses an initial belief/opinion 
modeled as a vector Xq G M*^. Agents are activated according to independent 
Poisson processes in continuous time.^ If agent a is activated at time t, her 
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Analogous versions of this model have been presented in the literature with agents' 
activations occurring in discrete time. 
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opinion jumps from its current value to a new value = u)X^_ 
where h is another agent sampled from A, and a; = 1 — (J G [0, 1] is a param- 
eter modeling how much agent a trusts the opinion of agent h. In general, 
the conditional distribution of h over the agent population may depend on 
the activated agent a (the support of such distribution representing the out- 
neighborhood of a in an underlying "social network" structure), while the 
parameter oj may depend on the interacting agents, a and 6, as well as on 
their current opinions, and X^_ . 

Fundamental theoretical issues concern the behavior of such models for 
large t and large n. Rather then in the single opinions' behavior, one is inter- 
ested in the emerging collective behavior of the population. Typical questions 
include whether a consensus is eventually achieved or rather disagreement 
persists, and, more in general, whether an asymptotic distribution of opin- 
ions exists, what it looks like, and how long it takes the system to approach it. 

The simplest case is when the Poisson processes are all of unitary rate, the 
conditional distribution of the observed agent is uniform over A whichever 
agent is activated, and the parameter uj is fixed and the same for all agents, 
independently of their current opinions. In this case, the model is linear and 
can be studied in full detail: it corresponds to the asymmetric gossip model 
in [18]. The basic fact is that (if u e]0, almost surely, all X" converge, 
as t — )• +00 (and for any fixed n), to a consensus random value ^ which has 
expected value E(^) = n~^Ylia-^(i- Convergence is exponentially fast [17]: 



E 



n-'y\xt-i{' 



a 



<2n~i^|X^|2exp(-Ct), 



where C = —nln(l — 2n ujuj — 2n-^uj^). The variance of ^ can be estimated 
as 



Var[e]<^^n-iy|Xol2 



Moreover, using the techniques in [18], one can easily prove a concentration 
result of type 

P(|Xf - E(Xf )| >e)< exp{-Ke'^n/t). 

Essentially, this shows that, as n grows large, and t/n tends to 0, each agent's 
opinion X^ concentrates around a deterministic dynamics converging to E(^) 
as exp(— 2a;ti;t). It is this type of results which we would like to extend to 
more general models. 

A particularly interesting setting is the homogeneous-population, state- 
dependent model, that is, when the parameter lo is independent of the iden- 
tity of the interacting agents, but does depend on their current opinions. 
The case 

(1) ^-a;(X,-,X,_)-|^^ if\X-_-X^_\>R, 
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where R> 0, and ujq G ]0, 1[, is known as the Deffuant-Weisbuch model [14, 
22, 23] of bounded confidence opinion dynamics: agents with opinions too 
far apart do not trust each other, hence they do not interact. Another case 
is the so-called Gaussian interaction kernel 

(2) uj = uiX^- , ) = uo exp(- - X^^ l^a^), 

a similar form of which was considered in [15]. Observe that, in these mod- 
els, the dynamics of the network and of the opinions become intertwined. In 
fact, these models are nonlinear and, to the best of the authors' knowledge, 
the only theoretical result [23] is that, if w E {0} U [coq, 1] for some > 0, 
each X^ converges, as t grows large, to a limit random value Numerical 
simulations show the asymptotic emergence of opinion clusters whose num- 
ber and structure depends on the initial condition but seems to be stable 
for large n. However, there is no theoretical result regarding concentration 
and scaling limits for any state-dependent model. 

On the other hand, for the case when the parameter co depends on the 
agents, as well as on their opinions, no theoretical result is available in the 
literature. Some of these heterogenous models have been considered in [19, 
25, 26, 31] where, though, only numerical simulations have been presented. 
Such heterogeneous population models are going to play a very important 
role in opinion dynamics because they are the natural model to represent 
more realistic populations with agents having different attitude to change 
opinion, and interacting only with agents in their social neighborhood. 

In this paper, we study general state-dependent gossip models for large n. 
We shall consider both the case of a homogeneous population, and of a het- 
erogeneous one consisting of two classes of agents: "standard" agents, which 
keep on updating their opinions as a result of interactions with the whole 
population, and "stubborn" agents whose opinions are never updated [1]. 
The latter case can be modeled as a homogeneous population model with 
an exogenous input describing the influence of the stubborn agents' opinions 
on the standard agents' ones, and interpreted as a, typically heterogeneous, 
"influential environment." We believe that many more general heterogeneous 
models can be studied with our approach. This will be done in a forthcom- 
ing paper where also models with interactions of nongossip type will be 
considered. 

In our analysis, we shall adopt an Eulerian viewpoint: instead of study- 
ing the evolution of the single agents' opinions, we shall neglect the agents' 
identities, and study the dynamics of the corresponding empirical opinion 
densities. We shall argue that the deterministic mean- held dynamics ob- 
tained in the limit of large n is governed by an ordinary differential equation 
(ODE) on the space of probability measures over the opinion set, presented 
in Section 2.2. As proved in Section 3, the initial value problem associated 
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to the mean-field dynamics always admits a unique global solution. More- 
over, at any finite time, its solution is absolutely continuous with respect 
to Lebesgue's measure, provided that so does the initial condition, and that 
some mild technical conditions are satisfied by the interaction kernel. 

The asymptotic behavior in time of the mean-field dynamics is analyzed 
in Section 4 for the state-independent heterogeneous case, and for the gen- 
erally state-dependent homogeneous case. In both cases, we prove weak con- 
vergence to an equilibrium distribution, which typically does not consist of 
a single Dirac's delta. For the state-independent heterogeneous model, we 
show that the equilibrium opinion distribution is independent of the initial 
condition, and is uniquely characterized by its moments, which can be com- 
puted by recursively solving a lower-triangular infinite linear system. On 
the other hand, we prove that the equilibrium opinion distribution in the 
bounded-confidence model is a convex combination of Dirac's deltas. Such 
deltas represent opinion clusters, and their number and position depend on 
the initial condition.^ Our results provide fundamental insight into two ba- 
sic mechanisms which have been proposed by social scientists in order to 
explain persistent disagreement in the society [4], namely heterogeneity of 
the social environment, and homophily leading to global fragmentation. 

Finally, in Section 5, we prove that the finite-population stochastic system 
concentrates around the deterministic mean-field dynamics, as the popula- 
tion size grows, at an exponential probability rate. We apply here a mar- 
tingale argument (see, e.g., [32] for the finite-dimensional case) and obtain 
a result in the Kantorovich-Wasserstein metric [3, 30]. The technical as- 
sumption in our results is that the, possibly stochastic, dependence of the 
weight (jj on the opinions is Lipschitz-continuous. Hence, the case (1) is not 
covered by our theory. This is not a relevant drawback since one can consider 
suitable Lipschitz approximations of (1); on the other hand, we believe that 
this is just a technical question and that the result should remain valid for 
a larger class of functions. 

We conclude this section with a brief overview of some related work. 
A special instance of the measure-valued ODE analyzed in the present pa- 
per has already been proposed in [5] for probability densities (in this case 
it becomes an integro-difFerential equation), but with no proof of either 
well-posedness or concentration of the stochastic finite system. In [6, 7, 10], 
deterministic, bounded-confidence, opinion dynamics models with possibly 
a continuum of agents have been studied both in discrete and continuous 
time. In particular, the continuous-time opinion dynamics studied in [10] is 
governed by a partial differential equation in the space of probability mea- 
sures, while the work [6] deals with the equivalent dynamics, in dimension 



^Proofs of similar results showing convergence of various variants of the bounded con- 
fidence opinion dynamics to opinion clusters have appeared in [6, 7, 10, 23, 25]. 
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one, of the cumulative distribution functions. In both works, the agents' 
opinions have continuous trajectories, and the corresponding generator of 
the opinion density dynamics is local. In contrast, in the model analyzed in 
the present paper, the opinion trajectories are discontinuous (in fact, piece- 
wise constant), and the induced mean-field dynamics is driven by a nonlocal 
operator. As shown in Section 4.1, the bounded-confidence mean-field dy- 
namics studied here has a qualitatively similar behavior to the solution of 
the partial differential equation of [6, 10]. It is also worth mentioning the 
work [21], where mean-field limits have been analyzed for the flocking dy- 
namics of Cucker and Smale [12, 13]. Finally, only at the end of their work 
the authors have become aware that an approach very similar to the one in 
this paper has been undertaken in [27], based on results in [20]. 

2. Problem setting and main results. In this section, we formally state 
the model and present our main results. 

Before proceeding, let us establish some notation to be followed through- 
out the paper. For x, y € W^, for some d € N, |x — y| and x ■ y will denote their 
Euclidean distance, and scalar product, respectively. The indicator function 
of a set A will be denoted by 1a, that is, lyi(x) = 1 if x G ^4, and tA{x) = if 
X ^ A. Given an open subset X C M*^, we denote by B{X) its Borel a-algebra, 
and hy J^{X) the space of finite signed Borel measures on X, equipped with 
the topology of weak-* convergence, while C M[X) denotes the 

closed convex cone of Borel nonnegative measures, and V{X) C Ai~^{X) the 
simplex of probability measures over X . The space of real- valued continuous 
bounded (resp., compact-supported, vanishing at infinity) functions on X, 
equipped with the supremum norm Hv^Hoo := sup{|(^(x)| : x € X}, will be de- 
noted by Ch{X) [resp., Cc{X), Cq{X)\. The Dirac delta measure centered in 
X ^ X will be denoted by 5x- For /i G M{X) and ip G Cfe(Af), we shall write 
ifJ-jf) for the integral / (p{x) dfi{x), with the convention that, whenever not 
explicitly indicated, the domain of integration is assumed to be the entire 
space X. The total variation of /i G M{X) will be denoted by ||/i||. The sym- 
bol A will denote Lebesgue's measure on X, <^ X will stand for absolute 
continuity, and dfi/dX for the Radon-Nikodym derivative, of fi with respect 
to A. Finally, we shall denote by Vi{X) := {fi G ^{X) : f \x\ dfi{x) < +00} 
the metric space of probability measures with finite first moment, equipped 
with the order- 1 Kantorovich-Wasserstein distance. The latter is defined by 
Wid^jiy) := mf{JJ \x — y\ dS^{x,y)}, where the infimization (which is in fact 
a minimization [3, 30]) runs over all couplings of fi and z^, that is, joint 
probability measures £ V{X x X) having marginals given by /x, and u, 
respectively. 

2.1. Stochastic models of continuous opinion dynamics. The present pa- 
per is concerned with continuous opinion dynamics systems. Agents belong 
to a finite population A of cardinality |^| = n. At time t G M"^ each agent 
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a £ A maintains an opinion G X, where X QMf^ is an open set. The 
vector of the opinions will be denoted by Xt := {Xf :aeA}e X-^. 

We shall assume the initial opinions Xq to be a collection of independent 
and identically distributed random variables, the law of each X^ given by 
some G ViW^). The trajectories of the opinion profile vector {Xt : t G M^} 
are right-continuous and evolve according to the following jump Markov 
process: Agents have clocks which tick at the times of independent rate-1 
Poisson processes. If her clock ticks at time t, agent a updates her opin- 
ion X^_ to a new value X^ which depends on the observation of the current 
opinion of some other agent and of her own one. In particular, she observes 
the opinion of some other agent h sampled uniformly from A, and then 
updates her opinion to a random value X", which has conditional proba- 
bility law k{-\X^_,X\_). Here k{-\-,-) is a stochastic kernel, that is, for all 
x,y G Af, is a probability measure on X, and i— )■ K{B\x^y) is 

a measurable map from X x X io [0, 1], for all measurable sets B <^ X. We 
shall refer to k as the interaction kernel of the model. We shall assume that 
the above stochastic process is defined on some filtrated probability space 
(r2,{J"t}jglR+,P), and denote by = Tq < Ti < T2 < • • • , the times at which 
some opinion update occurs (strict inequalities holding almost surely). Ob- 
serve that {Tfc+i — Tk'.kG Z"*"} is a family of independent rate-n Poisson 
random variables. 

In most of the models considered in the literature, the interaction ker- 
nel is a convex combination of type: K{-\x,y) = aK^{-\x,y) + aK^{-\x) where 
a = 1 — a G [0, 1] and where tv{-\x,y) is a probability measure concentrated 
on the interval connecting x and y while n^{-\x) is a probability measure 
concentrated on the segment connecting x to some random point z. More 
specifically, C M'^ is a convex open set containing the support of the initial 
condition, and there exists two scalar stochastic kernels ^*(-|-, •) and 0^{-\-, •) 
from X X X to [0,1] such that 



where = 1— uj,v = l — v and tp G V{X). This models a situation in which, 
with probability a, the activated agent updates her opinion towards a con- 
vex combination of her current opinion x and the opinion y of an observed 
agent. The weight u in such a convex combination measures the confidence 
that the activated agent has on the observed opinion of another agent, and 
is assumed to depend, through the stochastic kernel 0^{-\-, ■), on both the ac- 
tivated and the observed agent's opinions, x and y. On the other hand, with 
probability a, the activated agent observes an external signal z, sampled 
from a probability distribution tp, playing the role of an exogenous source 
of influence, or influential environment, and she updates her opinion toward 
a convex combination of her current opinion x and the observed signal z. 
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The dependence of the weight v of such convex combination is captured by 
the stochastic kernel 9^{-\-,-). A useful equivalent way to characterize the 
interaction kernel k described above is through its action on continuous test 
functions: 



for all ipeCoiX). 

Example 1 (Gossip model with heterogeneous influential environment). 
Assume that the stochastic kernel has the form (3), with constant 

weights: 



for some fixed confidence weights u},v £ [0,1]. This models a homogeneous 
population whose opinion dynamics alternates internal gossip updates to 
interactions with a static, external influential environment. Internal gossip 
steps occur with probability a, and involve a uniformly sampled agent a 
updating her opinion to a convex combination, with trust parameter w, of 
her current value and the one of another uniformly sampled agent b. Inter- 
actions with the external environment occur with probability a, and involve 
a uniformly sampled agent a updating her opinion to a convex combina- 
tion, with trust parameter w, of her current value and an external signal z 
sampled from a static distribution 'ip{dz). This model has been analyzed 
in [1] for finite, possibly inhomogeneous populations. The mean-field limit 
of this model, with homogeneous population, will be analyzed in detail in 
Section 4.1. 

Example 2 (Bounded confidence opinion dynamics). Consider the case 
when k{-\-,-) is in the form (3) with a = 1, and trust parameter distri- 
bution 6'^{-\x,y) supported on [0,a;o] for some ujq G [0,1[. The case when 
6^{-\x,y) = 6uj{x,y)^ where u}{x,y) is a nonincreasing function of the dis- 
tance |x — ?/| can be consider to model a homophily mechanism whereby 
agents are more likely to interact with others which have similar opinions. 
In particular, the case when uj{x,y) = for all \x — y\ > R, for some fi- 
nite R> 0, is usually referred to as bounded confidence opinion dynam- 
ics [14, 23], and the minimum such R as the confidence threshold. The spe- 
cial case uj{x,y) = u;o]1[o,_r](|2; — y\) corresponds to the Deffuant-Weisbuch 
model [5, 14, 24]. The mean-field limit of the bounded confidence opinion 
dynamics model will be analyzed in detail in Section 4.2. 



{K{-\x,y),(f) = a I ip{uox + Loy) d9^{Lo\x,y) 



(3) 




+ a 



9\-\x,y) = 6U-) 



r(.|x,y) = <5„(-) 
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While for most of the resuhs of our paper we shah not need the interaction 
kernel k to have the specific form (3), we shall focus on kernels of this 
form in Section 4 when proving asymptotic properties of the solution of the 
corresponding measure-valued ODE. 

Remark 1. The models considered in the cited literature often assume 
the interaction to be symmetric: when agent a is activated and interacts with 
agent b, both agents update their opinions. This symmetric model may be 
more suitable in certain applicative contexts, the asymmetric one in some 
others. However, while for finite population sizes some of the properties 
of the two models differ (e.g., in the symmetric model the average of the 
opinions is preserved, while this is not necessarily the case for the asymmetric 
model [18]), all the results and proofs of this paper hold, with minor changes, 
for the symmetric model too. 

2.2. The Eulerian viewpoint and main results. As the main interest is 
in the global behavior of the opinion dynamics system, rather than on that 
of the single agents' opinions, it proves convenient to adopt an Eulerian 
viewpoint, studying the evolution of the empirical densities of the agents' 
opinions. Formally, this is accomplished by considering the random flow of 
probability measures 

n ^-^ * 

This is a 'P(Af)-valued process whose trajectories are piecewise constant and 
right continuous. In particular, one has 

where {Mj. : k G Z+} is a P(A')-valued Markov chain. 

In order to describe the dynamics of the M^'s it is useful to consider the 
operator F : M+ {X) ^ M+ {X) , defined by 

F(/i)(i3):= j I K{B\x,y)d^x{x)d^Ji{y) ^BeB{X). 

Equivalently, one can write 

(4) {F{fi),ip) := j j j ip{z)dK{z\x,y)dfi{x)d^i{y) 

for all (f £ Co{X). When /i is a probability measure, then F(fi) may be inter- 
preted as the conditional distribution of the new opinion formed as a result 
of the first interaction occurring, given that the current empirical opinion 
density is fi. In fact, the opinions x and y of two agents a and b, ran- 
domly sampled, independently and uniformly, from the agent population. 
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have conditional joint distribution dfi{x) dfi{y) , and hence the new opin- 
ion z formed as a result of their interaction has conditional distribution 
dK{z\x, y) dii{x) d^{y). 

It is immediate to verify that 

E[(Mfc+i, (^)lMfc] = (1 - n-i)(Mfc, + n-i(F(A4), V9) 

for all if ^Cq{X) and k G One may rewrite this in the form 

(5) (Mfc+i,(^) - (Mfc,^) =n-i(F(Mfc) -Mfc,(^) +n-i(Afc+i,(/p), 
where the random signed measure A^+i satisfies 

(6) E[Afe+i|JTj=0, ||Afc+i||<n||Mfc+i-Mfc|| + ||F(il/4)-Mfc||<4. 

Equation (6) implies that {(A^, 93) : G N} is a sequence of bounded martin- 
gale differences, which can be thought as "noise." This suggests to think of 
the stochastic process {M^ : k G Z+} as of a noisy discretization, or Euler ap- 
proximation in the numerical analysis language, of the probability-measure- 
valued ODE 

(7) ^Att = F{iit) - fJ-t 

with stepsize 1/n. We shall refer to a solution of (7) as the mean- field dy- 
namics of the system. 

More precisely, we shall define a solution of (7) to be a family of probability 
measures {fit :t £[0, +00)} such that, for every function 99 G Co(Af ), the real- 
valued map 1 1— 7- (/ii , ip) is differentiable on M"*" , and satisfies 

(8) '^^ " (Fif^t),^) - {fJ-t, V>) 

for every t > 0. The main result of this paper, stated below, guarantees 
that (7) admits a unique solution {/it}, and that the stochastic process {/i"} 
concentrates around {nt} exponentially fast in n. 

Theorem 1. Let ii£V{X) he arbitrary. Then: 

(a) There exists a unique solution {fit -t G M"*"} of (7) with initial condi- 
tion fio = fi; 

(b) If X (^W^ is bounded, and the stochastic kernel k is globally Lipschitz 
continuous as a map from X x X to Vi{X), then, for every r G (0, +oo), for 
sufficiently small e > and sufficiently large n G N, it holds 

F{sup{Wi{fi'l,fit):te[0,T]}>e)<exp{-K£\), 

where K is a positive constant depending on X , k and r only. 

Points (a) of Theorem 1 will be proved in Section 3.1, while point (b) will 
be proved in Section 5. Additional properties of the solution of the initial 
value problem associated to (7) will be studied in Section 3.2, while Section 4 
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will present an analysis of the behavior of the mean-field dynamics for the 
model with heterogeneous influential environment, and for the bounded- 
confidence opinion dynamics. 



3. Well-posedness of the measure-valued ODE. In this section, we shall 
first prove point (a) of Theorem 1, that is, that the initial value problem 
associated to the ODE (7) admits a unique solution. Then, under further 
technical assumptions, we shall show that, if the initial measure //q admits 
a density, so does the solution fit at any finite time t. 



3.1. Weak solutions. To start with, we extend the ODE to the space of 
signed measures In order to do this, we need to extend the opera- 

tor F and introduce another operator G in the following way. For // G A4{X), 
put 

(9) Fif,):=Fi^L+), G{^,):=^,+ {X)f,, 

where n = jj,'^ — denotes the Hahn- Jordan decomposition of ;U G Ai{X). 
It is not hard to check that both F and G are locally Lipschitz continuous 
with respect to the total variation norm, that is, for every bounded set 
O C A^(A'), there exist nonnegative constants Kf,Kg such that 

||F(/ii)-F(/X2)|| <i^F||/il-/U2||, 

(10) 

||G(/ii)-G(/Z2)||<i^G||/Xl-/X2|| 

for all Hi,H2 G ©• Moreover, 

(11) F(/.)(A') = G{f^){X) = ,x{Xf V;u G M^{X). 

In the following, we want to study the well-posedness of initial value 
problems associated to the measure-valued ODE 

(12) J^^xt = F{^lt)-G{^lt), 

where (12) means that, for every G Cq{X), the real-valued map 1 1— )• (//t, (p) 
is differentiable on M+, and satisfies ^{fit,f) = {F{^it),'p) — {G{iJ,t),(p) , for 
every t > 0. We shall refer to such a {nt : i > 0} as a weak solution of (12). 



Proposition 1. Suppose that F,G:M{X) ^ M~^{X) satisfy proper- 
ties (10), and (11)- Then, for every ^^M.'^{X), there exists a unique so- 
lution {nt : t G M"*"} C 7W+(Af) to (12) such that fj,o = fx. Moreover, fJ.t{X) = 
fi{X) for every t>0. 



ds 
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Proof. For r G (0,+oo), let C([0,r], be the space of continu- 
ous curves in M^{X) equipped with the sup norm ||{/Uf}||T- := sup{||/it|| :t G 
[0,r]}. Given a curve {fj,s} G C{[0,t],M.{X)), and a bounded measurable 
function &Co{X), define 

m{fis})t,f):={fi,v)+ [ {F{fis),v)ds- [ {G{fis),v) 

(13) ^° 

Vt G [0,t]. 

Observe that (12) with the initial condition /xq = /U is equivalent to 

(14) {fit,^) = {H{fis})u^) Vv9GCo(A'), t>0. 

Notice that, for every t G [0, r], ^{{ns})t can be seen as the difference of two 
bounded linear positive functional on Co(Af), so that <&({^s})j G 
Moreover, the map $({^s})j is continuous over [0,r], since 

\M{f,s})t+e-HM)t\\= r'\\G{fis)\\ds+ \\F{fis)\\ds 

(15) 

<e[\\{G{f,s)}\\r + \\{F{f^s)}U 

Therefore, the operator ^ takes values in C{[0,t],M{X)). Now, let us con- 
sider := {z^ G : ||i^|| < 2||/i||}, let Kf,Kg be the Lipschitz constants 
relative to © of F, and G, respectively. For every G 0, (11) and (10), imply 
that 

(16) ||F(z/)|| < ||F(z/) - F(^)|| + ||F(/i)|| < 4KfM. 
Similarly, 

(17) \\G{i^)\\<AKgM\. 

Define now the set S := {{fj-t} G C{[0,t], M{X)) : =fJ-,Ht& 6, Vt G [0,r]}. 
For all {fJ,t} G S, using (16) and (17), and arguing like in (15), we obtain 

(18) ||«5({M)llr<(l+4ri^)||/i||, 

where K := Kp + Kq- Moreover, if both and {I't} belong to S, then, 

IMM) - H{^t})\\r = sup [\\\F{^Ls) - F{us)\\ + \\G{us) - G{^^s)\\) ds 
(19) 

<rK||{M-{^ailr. 

We now assume to have chosen r G]0, ^^J. Then, by (18), <^(5) C S and, 
by (19), $ is a contraction of S. Hence, by Banach's fixed point theorem 
there exists a unique fixed point of <l> in 5. As observed, such a fixed point 
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corresponds to a solution {fJ-t} of the ODE (12) for t G [0,r], with the ini- 
tial condition fiQ = fi. We now show that indeed fit G for t £ [0,r]. 
By contradiction, assume that there exists B G B{X) such that /it(-B) < 
for some t £ [0,r], and let t* := sup{s € [0,t]: iis{B) > 0}. By continuity, 
fit*{B) = while fJ-s{B) < for all s £]t*,t]. This implies that 



Fifis)iB) - GifisKB) > -fitWfis{B) > Vs G ]t*,t]. 



MB)= / iFifi,){B)-GifMs)iB))ds>0, 



But then 



which is a contradiction. Hence, fit G M'^^X) for t E [0,r]. Notice moreover 
that, because of property (11), fJ-ti'^) = fJ'i^) for all t G [0,r]. Finally, a stan- 
dard induction argument allows one to extend the existence and uniqueness 
of the solution to the whole interval [0, +oo). □ 

Notice that, when considering an initial condition ^uq € V{X), the solution 
of (12) satisfies /Uj G 'P(Af) for all t, thus proving point (a) of Theorem 1. 

3.2. Probability density solutions. We shall now investigate on the exis- 
tence of density solutions when the initial condition /iq is absolutely contin- 
uous with respect to Lebesgue's measure. 

Given the interaction kernel k(-|-,-), and a nonnegative measure fi in 
M+{X), we put 



i^i{lJ-){B\y) ■■= / K,{B\x,y)dn{x), 



K2{n)iB\x) := / K,{B\x,y)dn{y) 



(20) 



for all B G B{X), x,y £ X. The following result characterizes regularity prop- 
erties of the solution of the initial value problem associated to the ODE (12). 

Proposition 2. Assume that /io <^ A, and that 

(21) ^<A =^ Ki{n){-\y),K2{n){-\x) <^X yx€X,yyeX. 

Then, fj.t <^ X, for all t G [0, +oo) . Moreover, if there exists C G (0, +oo) such 
that, for all fi<^ X, 



(22) 



dK2{fl){-\x) 



dX 



<C 



dfi 



dX 



yx€X, 



then, the density ft = dfj,t/dX satisfies the estimation 



(23) 



IIAIIoo<||/o| 



Vt G [0,+oo). 
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Proof. For every finite time t G [0, +00), consider Lebesgue's decom- 
position fj,t = fj.f + nf , where <C A, and fif and A are singular. It follows 
from (21) that, K2{fJ.t){-\x) < A for all x€X. Then, for any B £ B{X) such 
that \{B) = 0, one has 



+ 



K{B\x,y)dn^{x)dfit{y) = J dK,2{l4){B\x) d^it{x) = {). 
Similarly, one can show that JJ K{B\x,y) d^Kx) d^Ky) = 0. Hence, 
F(/it)(5) = j j K{B\x,y)dnt{x) dfitiy) 

K{B\x,y)dfj.^{x) dnt{y) + J j K{B\x,y) dnt{x)dn^{y) 

K{B\x,y) dfif{x) dfit{y) 

K{B\x,y)dfj.t{x)dfit{y) 

for all B G B{X) such that A(-B) = 0. This readily implies that /i| satisfies 

Since /Uq = by assumption, it follows that = for all t>0. 

Assume now that (22) holds true. For any if G Cc{X), Holder's inequality, 
and (22) imply that 



ip{z) dK{z\x, y) d^tix) dfit{y) 



^(z)^^^^;^M^dA(z)d/.,(x) 



< 



dX 

dK2ifit){z\x 



dX 

< Cll/tlloollv'lll- 

It follows that, for all nonnegative- valued ip £Cc{X) 



ifWi dnt{x) 



ip{x)ft{x)dX{x) = J (p{x)fQ{x)dX{x) + J {{F{fis),'^) - {tJ's,'^))ds 
< ll/olloolblli + / {F{^J's),^p)ds 



< 



m\i \ ii/o 



+c j\\fs\\^dsy 
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Then, by the isometry of L^{X) with the dual of L^{X), the fact that ft is 
nonnegative valued, and the density of Cc(Af) in L^{X), one gets that 

ll/tlloo = supjy (p{x)ft{x)dx:(f e L^{X), \\(f\\i<l^ 

= supj j (p{x)ft{x)dx:(f e Cc{X),(f> 0, II^jIIi < l| 



< \\fo{x)\\oo + C / WfsWoods. 
Jo 

By Gronwall's lemma, this readily implies (23). □ 



The technical condition on the stochastic kernel k is actually verified in 
many important cases encompassing the bounded confidence dynamics (1) 
as well as the Gaussian interaction model (2). 



Corollary 1. Assume that the interaction kernel k is the form (3) with 
^'{■\x,y) = 5^i\x-y\) and 9''{-\x,z) = where a;:M+ [0,tJo], ^^0 G 

[0,1[, and f:M+ — > [0, t^o], ^^o £ [Ojll; (^^^ both nonincreasing and piece- 
wise . If fio<^ X, then fif^ X, for allt G [0, +oo) and the relative densities 
satisfy condition (23). 

Proof. We shall show that the conditions of Proposition 2 are satisfied 
in this case. Fix y ^ X and consider the function x i— t- a)(|x — y\)x + oj{\x — 
y\)y- The assumption on u ensures that it is an invertible transformation in x 
and a simple geometric consideration shows that the inverse has the form 

x = g{w,y) = y + a{\w-y\){w-y), 

where a:M'^ — )• M"*" is such that a{Co{t)t)Co{t) = 1 for all t>0. The function 
r](t) = u){t)t is strictly increasing, hence invertible and we can thus write 
a{s) = [u){ri^'^{s))\^'^ . a is thus also a piecewise function as well as g{-,y) 
whose Jacobian can easily be shown to be 

Dwg{w,y) = a{\w - y\)I + Va{\w - y\){w - y)K 

Straightforward computation show that D^g is bounded in the pair {w,y). 
Similarly, the function x i— )• v{\x — z\)x -\-uj{\x — z\)z admits an inverse in x, 
X = h{w,z) whose Jacobian D^h is bounded in the pair {w,z). 

Then, if fi is absolutely continuous with density /, one has for all non- 
negative real- valued 99 G L^{X), and y ^ X, 



{i^i{lJ'){-\y),V>) =ct J f{u}{\x-y\)x + uj{\x-y\)y)f{x)dX{x) 

+ a ip{v{\x — z\)x + v{\x — z\)z)f{x)dX{x)d'ip{z) 
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= aj ip{w)f{g{w,y))\D^g{w,y)\dX{w) 

+ a J J ip{w)f{h{w,z))\Dy^h{w,z)\dX{w)dip{z) 

<Ci\\f\\M\i, 

where Ci := a\\Dwg{w,y)\\oo + a\\Dii,h{'w,z)\\oo- Similarly, one shows that 
there exists some constant C2 > such that {K2{lJ-){-\x),ip) < C2||/||oo||v'lli) 
for all nonnegative valued E Li{X) and x G ^. As a consequence, i^i{-\y) ^ 
A, for all y S A' and i = 1,2, and (22) holds. Therefore, the claim follows from 
Proposition 2. □ 

4. Behavior of the mean-field dynamics. This section is devoted to a dee- 
per analysis of the ODE (12) for the state-independent gossip model with 
heterogeneous influential environment, and the bounded-confidence opinion 
dynamics, respectively. In particular, we shall investigate the limit behavior 
as t grows large, showing that, in both models, fit converges weakly to an 
asymptotic opinion measure. The behavior of the two models, and their anal- 
ysis, however, differ substantially. For the state-independent gossip model 
with heterogeneous influential environment, the ODE governing the mean- 
field dynamics is linear, and can be analyzed by iteratively solving the lower- 
triangular linear system of ODEs governing the various moments behavior. 
In this case, the asymptotic opinion measure is independent of the initial 
value, it is characterized by its moments, and is absolutely continuous if so is 
the influential environment. In fact, one could show that the corresponding 
finite population Markov process is ergodic. In contrast, the ODE govern- 
ing the mean-field dynamics of the bounded confidence model is nonlinear, 
and convergence is shown by a Lyapunov argument. The asymptotic opinion 
measure is given by a convex combination of deltas, whose number and po- 
sition typically depends on the initial condition. Indeed, the corresponding 
finite population Markov process is typically not ergodic, in this case. 

4.1. Gossip model with heterogeneous influential environment. We start 
by analyzing the case when the stochastic kernel k{-\-,-) has the form (3), 
with constant weights: 9^{-\x,y) = 5i^{-), 9'^{-\x,y) = 6y{-), for some fixed 
uj,v € [0, 1]. Throughout this section, we shall assume an exponential bound 
on the moments of both ^lq and tl^, that is, 

(24) sup( / \x\'' dfioix) ] <+oo, sup( / \xfd'il;{x)] <+oo. 

Clearly, (24) is automatically satisfied when X is bounded. Let us fix some 
z G M'^, and consider the z-weighted moments of fit and ip, respectively, 

^^'^ - f{x-z)''dfit{x), nf^ := l\x-yfdij{y), keZ+. 
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The following result characterizes their evolution in time. 

Proposition 3. The z-weighted moments satisfy 
(25) ^mj^^ = au(n(^) -mf)), 



(26) ^mf^ = -7femf^ + /fc(mf\...,mf-^^) + at;V=), k>2, 
where 



/fc(mj^\...,mf"^^) := V ( ] {au,^uj''-^my'm\''~" + av^v^-^mf'n 



Proof. For the first moment, one has 

d_ 
~dt 



^m[^^ = a I I {{ujx + uiy) ■ z) djitix) dixt{y) 



(1) 



^ ^ j j (('^^ ~^ "^y) ' ^) dfJ-ti^o) dip{y) — 

= avn^^^ — avm[^\ 
which proves (25). For k>2, one has 

(Ox ■ z + ujy- z)^ dfxtix) dfitiv) 
= (w^ + a;^) j\x ■ z)^ d^it{x) 



and, similarly, 



{{vx + vy) ■ z)^ djJtix) dipiy) 



i=i ' 
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From the two identities above, it follows that 

d (k) 
dt 

fe , / N, , , / N (fc) 



ml ' = a I I {{u}x + uy) ■ z) ' dfitix) d^tiv) 



+ ay j {{vx + vy) ■ z) dfit{x) dil){y) — ml 

which proves (26). □ 

Example 3. In the special case when a = 1, namely when there is no 
influential environment, we obtain from (25) that ^ J xdpLt{x) = 0, so that 
the first moment is constant. On the other hand, the variance 



vt :-- 

satisfies ■^vt = —2uju)Vt. Hence 



X - I ydi^Loiy) 



dfitix) 



that is, fit converges to a delta centered in the average initial opinion expo- 
nentially fast in t. 

We now focus on the limit as t — )■ +oo for the general case. An inductive 
argument proves the following result. 

Lemma 1. Assume a < 1. Then, for every z G W^, the z-weighted mo- 
ments of Ht satisfy 

(27) lim mP=m^^\ k£Z+, 

t— >oo 

where m!'^ can he recursively evaluated by 

(28) mW := n^') , ^(^+1) = 7,^^i[/fc+i(mW, . ..,mg))+ av'n(>^+% 
Proof. For k = l, the solution of the ODE (25) is easily found to be 



(29) ml'> = e-""*m^'^ + (1 - e 

so that equation (27) clearly holds. Moreover, assume that equation (27) 
holds for every A; G {1, . . . ,j — 1}, and define ■= fj{m[^\ ■ ■ ■ ,mj"' ^^) for 



t € [0, +oo]. Then, the continuity of fj implies that limt_s.oo xF"* = Xoo^- Solv 
ing the ODE (26) gives 

(30) mp'^ = / e-^^ (xF + av^n^^^) ds + e'^^^m!'^'^ . 
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Clearly, the second addend of the right-hand side of (30) converges to zero 
for t — )• oo. On the other hand, the convergence of Xt imphes that 

lim [ e-^^^*~'\xt^ + av^n^^'^)ds = {x^£ +av^n^^^) lim [ e~^^^^~'Us 

The foregoing, together with (30), implies the claim. □ 

We are now in a position to prove the following result for the convergence 
of f-H- 

Proposition 4. Assume that (24) holds. Then 

lim Ht = fJ'oo, 

t— >oo 

weakly, where fi^o & 'P{X) is uniquely characterized by its moments . 

Proof. It follows from (24) that there exists some finite M £ M"'' such 
that 

(31) |m[,*^^| < |z|*^M^ \n^^^\<\z\''M'' 

for all 2 G M'^ and A; G N. Now, an inductive argument shows that 

(32) |?nf ^1 < I^I'^M'^ VtG [0,+oo], zGM'^ 

for all A; S N. In fact, (29) and (31) immediately imply that (32) holds for 
k = l. Moreover, if (32) holds for all /c e {1, . . . , j - 1}, then (30) and (32) 
give 

|mp)| < / e^^^^^^'\\fj{m[^\...,rn'f~^^)\+av^\n'^^^\)ds + e-"'^^\m^^^\ 



< [ e^^^^^-'^^M^zP-fids + e-^^^M^lzP 



7j 







Let us consider the characteristic functions (j)t{z) := f exp{iz ■ x)dnt{x) 
and, for k G Z+, define at{k) := i^m"^^ /k\, b{k) := M^\z\^ /k\, and observe 
that X]fc6Z+ ^(^) ~ exp(M|z|). One has that 

<t>t{z)= j J2 ^'^kf dfit{x)= ^ f {x ■ zf dfitix) = 
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where the exchange between the series and the integral is justified by Lebes- 
gue's dominated convergence theorem, since 



0<k<n 



0<k<n 



Moreover, observe that, since |of(A;)| < b{k), another apphcation of Lebes- 
gue's dominated convergence theorem gives 

lim (ptiz) = lim at{k) = aoo(^) =: (pcx^iz). 



fcez+ 



fcez+ 



Hence, (ptiz) converges pointwise to (pooiz), which in turn can be easily 
verified to be continuous at 0. Then, the claim follows from Levy's continuity 
theorem (fSl, Theorem 2.5.1). □ 



Observe that, for all a € (0,1), the limit measure fXoo is independent of 
the initial condition /Uq, and depends only on the influential environment ip, 
as well as on the parameters a, cj and v. Notice that the first moment sat- 
isfies = n^^h In contrast, if V' 7^ ^xq, it easily seen that 7^ n^^^ for 
/c > 2, so that in particular fi^o 7^ ip- On the other hand, it follows from (28) 
that, if V' 7^ ^xo) then the variance of /ioo is positive, so that fioo / ^xo - This 
result may be interpreted as showing that the presence of an heterogeneous 
influential environment prevents the population from achieving an asymp- 
totic opinion agreement. In fact, as shown in the following proposition, the 
asymptotic opinion distribution ^00 is absolutely continuous whenever so is 
the influential environment ip. 

Proposition 5. Assume ip <^ X. Then ^oo ^ A for all a G [0, 1). 

Proof. For ViX), 7 G [0, 1], deflne 7 := 1 - 7, and 

L^/ifi, v) G ViX), {Lyifi, v),Lp) ipi^x + 7y) (i/i(x) dviy) 

for every Lp G Cb(Af). Since is a rescaled convolution operator, and since 
tp <^ X, one has that Lyifi,7p) <^ X. Similarly, ^a;(/-fooi /^*oo) = Q;L(/i^, /i^), 
where is the singular part of ^00 • Combining this with the fact that the 
asymptotic measure satisfies 

Hoo = Finoo) = aL^ifi, fi) + aLyifi,'ip), 

one gets that 

M^(Af) = aiLUliL,l^L))iX) = aifi^iX))^. 
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t=2 



t=5 



t=10 



t=+cx 



dy/dx 



Fig. 1. Behavior in time of the ODE solution in d= 1, with initial condition /lo uniform 
over (0,10), heterogeneous environment dip{x) = exp{ — {l — {x — 3)'^)~^)l(^2,4,){x) dx, and 
parameters a = 0.5, lu = 0.5 and v — 0.5. The Radon-Nikodym derivates of the asymptotic 
measure /ioo, and of the influential environment (dashed) are plotted as a reference. 



Therefore, 

/i^W(l-a/i^(Af)) = 0. 
Since fi^{X) < 1 and a < 1, this necessarily impHes that /i^(Af) = 0. □ 

Figure 1 reports numerical simulations of the mean-field dynamics, when 
started from a uniform distribution over an interval, and influenced by an 
absolutely continuous environment. Coherently with Proposition 2, the so- 
lution remains absolutely continuous during its evolution. As t grows large, 
fit converges to a limit measure whose first moment coincides with that 
of ip, and which is absolutely continuous, as predicted by Propositions 4, 
and 5, respectively. Such a limit density may be interpreted as resulting 
from a tension between the aggregating forces represented by the first ad- 
dend in the right-hand side of (3), and the environment's influence captured 
by the second addend in the right-hand side of (3). 

4.2. Bounded confidence opinion dynamics. We analyze now the case 
when •) is in the form (3) with a = 1, and weight distribution 0{-\x, y) : = 
supported on [0,u;o] for some ujq G [0, 1[, and satisfying the symme- 
try assumption 

(33) e{-\x,y) = e{-\y,x) 

for all x,y S A'. The following result states weak convergence of fif 

Proposition 6. Assume that J |xp dfio{x) < oo. Then, there exists //qo £ 
V{X) such that 

lim fit = IJ'oo, 

t— >-oo 

weakly. 
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Proof. We start by proving that the second moment mf^ := J\x\'^ d^t{x) 
is a Lyapunov function for the system. Observe that, for ah x,y G M'^, 
we [0, 1] , a) = 1 — w, one has 

\x + uj{y - x)\'^ + \y + uj{x - = {uP' + u?){\x\'^ + \y^\) + AujCox ■ y, 

so that 

2ijjijj\x — = 2ti;;D(|j;p + |yp — 2x ■ y) 

= (1 - - + - 4u;wx • y 

= |rEp + |y^| — |x + a;(y — x)p — |y + ijj{x — 
From the foregoing, and the symmetry of 0{-\x,y), it fohows that 

^m^^'' = y jxp dF(//t)(x) — m|^^ 



(34) 



1 1 j j{\x + u:{y-x)\' 

+ |y + - - - {yl"^) d9{uj\x,y) dfitix) dfitiy) 

U}{1 -uj)\x- (i0(u;|x, y) d;U((x) d^tiu) 
<-(l-L^o)Tt, 



where 




Tt := uj\x-y\ de{uj\x,y) dfit{x) dfitiv)- 



Hence, in particular, -^ml < 0, so that is nonincr easing, and therefore 

(2') (2) 

convergent. Define := limt^oo nT-t observe that (34) imphes that 

rt 1 d 
hm / T^ds< lim / —m^?^ds 



t^co Jq t-5>oo 1 — LijQ Jq ds 



(35) = Hm 



(2) (2) 

uIq — ml 



t^oo 1 — Wo 
(2) (2) 

1 - Wo 

Now, for any smooth and compact-supported test function ip G C^(M'^), we 
can write 

(36) ip{x + uj{y - x)) - ip{x) = uj{y - x) • Vip{x) + r{x, y), 
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with |r(x,y)| < uj'^\y — xp<l> where $ := Moreover, again from the 

symmetry of 6{-\x,y), one has 

uj{y - x) ■ V(p{x) d9{u}\x, y) dfit{x) dfitiy) 

uj{y - x) • (V(/?(x) - Vip{y)) de{uj\x, y) dfitix) dfit{y) 
uj\y - x\\Vip{x) - Vv3(y)| de{u}\x,y) dfit{x) d^it{y) 
uj\x - d9{uj\x, y) dnt{x) dfitiu) 



(37) 



_ 1 
~ 2 
1 

< - 
- 2 



< — 
- 2 



< -Tt. 

- 2 

From (36) and (37), it follows that 



\{Fifit)-^lt,v>)\ 




< 




{ip{x + u}{y - x)) - ip{x)) d9{uj\x, y) d^t{x) dfit{y) 



u}{y - x) ■ V(p{x) d9{uj\x, y) dfitix) dfitiv) 



3$^ 
< — Tt, 
-2 




u}^\x - y\'^ d9{u}\x,y) dfit{x) dfj,t{y) 



so that 



/■* 3$ /■* 

lim / |(F(^s) - /is,<^)| ds < — lim / Tgds 



< 



3$ 



2(1 -wo) 



(m[,^^ -m(^)). 



Therefore, in particular, the limit 



lim {fit, f) = lim / {F{fis) - f^s,f) ds 



exists and is finite. From the arbitrariness of G C^(Af), it follows that 
converges in the sense of distributions. Finally, notice that, since the second 
moment is bounded, the family {/if :t G M+j is tight, hence /if converges in 
V{X). □ 

If we make the further assumption that the weight a; ~ 9{-\x,y) is almost 
surely strictly positive in a neighborhood of the diagonal {(x, x) : x G X}, we 
have the following characterization of the equilibrium points. 
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Proposition 7. Let R> be such that 

S{R) := ini{uj :supp{9{-\x,y)) C [oj, 1] \fx,ye X, \x - y\ < R} > 0. 

Then /ioo is a convex combination of Dirac's deltas centered in points sepa- 
rated by a distance not smaller than R. 

Proof. Assume by contradiction that there exist x*,y* G supp(yUoo) 
such that |x* — y*| < R. Then, one can find suitable neighborhoods A and B 
of X* and y* , respectively, such that \x — y\ < R for all x G A and y £ B. 
Hence, supp(6'(-|x, y)) C [5{R), 1] for all x G A and y £ B. Then, 



This clearly contradicts (35). □ 

It is worth stating the following simple, though important, consequence 
of Proposition 7, which, in particular, applies to the Gaussian interaction 
kernel (2). 

Corollary 2. Suppose that 



Then, /ioo = ^xo where xo = J xdiJ,o{x). 

Figure 2 reports numerical simulations of the mean-field ODE associated 
to the bounded-confidence model of Deffuant-Weisbuch, in dimension d= 1, 
starting from an initial condition uniform over the open interval (0, 10). Ob- 
serve that, as predicted by Proposition 2, the solution remains absolutely 
continuous, with bounded density, at any finite time t. It is possible to appre- 
ciate the effect of local aggregation forces, which first lead to the formation 
of two peaks around the opinion points x = 1,9, then of other two smaller 
peaks around the points x = 3,7, and finally of a smaller peak in x = 5. 
As t grows large, the opinion density converges to a convex combination of 
Dirac's deltas, as predicted by Proposition 6, separated by an inter-cluster 
distance of at least 1, as predicted by Proposition 7. A schematic represen- 
tation of the asymptotic opinion distribution, as studied in [5], is plotted as 
well, presenting some minor clusters between the major ones. The reader is 
referred to [5, 24] for extensive simulations of this model, and bifurcation 
studies for the asymptotic distribution. These results may be interpreted 





U {(x, y) : supp(e(-|x, y) C [uq, 1])} = A- x A". 



uio>0 
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t=o t=1 t=2 t=3 t=4 t=+~ 



Fig. 2. Behavior in time of the ODE solution in d= 1, with initial condition fj,o uniform 
over (0,10), and y) = (5i/2l[o,i] (|i — J/|) + 5ol(i,+oo)(|a; — 2/|) • A schematic represen- 

tation of the asymptotic distribution (t = +oo ), as presented in [5], is reported as well. 

as explaining how locally aggregating interactions modeling homophily can 
generate global fragmentation. 

We conclude this section by observing that arguments along the lines of 
the proofs of Propositions 6 and 7, combined with a standard martingale 
convergence theorem, can be used, for every finite population size n, to prove 
almost sure convergence of the stochastic system /i" to a random asymptotic 
measure /xj^, consisting of a convex combination of Dirac's deltas separated 
by a distance of at least sup{i? > : 6{R) > 0}. 

5. Concentration around the mean-field dynamics. In this section, we 
finally show that, as the population size n grows, the stochastic process {/x"} 
concentrates around the solution {fit} of the ODE (7), at an exponential 
probability rate. Throughout this section, we shall assume that Af C M"^ is 
bounded, with A denoting its diameter, and that the stochastic kernel •) 
is globally Lipschitz in the Kantorovich-Wasserstein metric, that is, that 

(38) Wi{K{-\x,y),K{-\x' ,y')) <^\{x,y) - {x',y')\ yx,x',y,y'£X 

holds for some finite positive constant Lp. Our first step consists in showing 
that the operator F inherits the Lipschitz property from the stochastic kernel 
k{-\-, •). The proof of the next result relies on the duality formula ([3], (7.1.2)) 

(39) Wiifi, v) = sup{(^, V?) - (i/, V?) : V9 G Lipi(A:')}, 
where Lip^(A') denotes the set of 1-Lipschitz functions on X. 

Lemma 2. If (38) holds, then 

Wi {F{fi) ,F{u))< LpWi ifi, v) V/i, u^Vi {X) . 
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Proof. First, observe that, for arbitrary function 93 G Lip^(Af), and 

Lp{z)dK{z\x,y) - j ip{z)dK{z\x',y') < Wi{K{-\x,y), k{-\x' ,y')) 

<^\(x,y)-{x',y')\ 

< -^(k-a; \ + \y-y |), 

by (39) and (38). For £ V{X), let ^ G Pi(A' x X) be their optimal cou- 
pling, that is, the one such that JJ \x — y\d^{x,y) = Wid^ju). Then 

V(.)<i.(*, 

ip{z) dn{z\x' ,y') dv{x') dv{y') 
(p{z){dK{z\x, y) - dK{z\x\ y')) d^(x, y) di{x\ y) 



^-^111 1 i\x-^'\ + \y-y'\)d^i^,y)d^{x',y') 



2 

= LFWi{fi,v). 

Hence, the claim follows by applying the duality formula (39) once more. 
□ 

Observe that there are three sources of randomness in the system: the 
empirical measure of the initial opinions //q , the update times {T^}, and the 
agents' interaction. The first two can be easily dealt with by appealing to 
the following classical large deviations results. 

Lemma 3. For all G V{X), e > 0, it holds that 

limsup?i^MogP(W^i(/i[J,/io) > e) < -e^ /2. 

n 

Proof. Sanov's theorem ([29], Theorem 2.14) and the Csiszar-Kullback- 
Pinsker inequality ([30], page 580) imply that 

liminf-ilogP(t^i(/iS,/io)>e) 

n n 

> \ni{H{v II iiq) : u G V{X),Wi{v, fio) > e} 
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> inf{i||i/ - /iof : V G W^i(i/,//o) > e} 

>eV(2A2), 

where iJ(z^ || [i) denoted the relative entropy, and the last inequality follows 
from the estimate W\{y^[x) < A||i/ — //|| ([30], Theorem 6.15). □ 

Lemma 4. For t G M+, let (;{t) := sup{k G Z+ : Tj, < t}. For all r € 1R+, 
a>l it holds 

limsupn"^ logP(sup{i - r<-(t) : < f < r} > e) < -eV^^^ 

n 

limsupji^^ logP(<j(T) > am) < — {a — l)^r. 

n 

Proof. The first statement follows, for example, from [29] , Theorem 5.1. 
The second one, for example, from [29], Example 1.13. □ 

We are now left with the third source of randomness, originated by the 
selection of the interacting agents, and their actual interaction. Observe 
that, in the right-hand side of the duality formula (39), one may restrict the 
supremization to the test functions (p belonging to 

Lipf :={^GLipi(3^):|^(x)|<A/2}, 

where 3^ is a hypercube of edge-length A containing X, and /i, are naturally 
identified as elements of V{y). The following result shows that the set Lip^ 
can be approximated in the infinity norm by not-too-large a set of functions. 

Lemma 5. Let CW^ be compact and convex. Then, for all e S ]0, A/2], 
there exists a finite set Tie C Lipf such that \ne\ < 2v|+i |3({A/£)(v^+i))'*^ 
and 

mm{\\h-ip\\:h£ne}<e Vv? G Lipf . 

Proof. With no loss of generality, we shall restrict to the case X '^y = 
[0,A]'^. We introduce a discretization operator Lipf —)• Lipf as follows. 
Let T] := e/{Vd+l/2) and define J := {0,1, . . . , [A/rj\}. For any ip G Lipf 
and j G J"^, let k{j) = i £ J iff ip{jr]) G [-1/2 + r]i, -1/2 + r]{i + Observe 
that, since ip is 1-Lipschitz, one has 

(40) Yl \ji-j'i\<^ =^ m-k{^)\<i. 

l<l<d 
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Then, define ^{ip) = h, by putting, for all x G ni<«<db'/^; Ui + 

h{x)= H f(k{i+6i)-kmxi-jiv)+vm-l + l)- 
i<l<d^ ^ 

Thanks to (40), one has that G Lip^ for all G Lip^- Moreover, for 
aU j G J'^, one has |$((^)(jr/) - (/?(j??)| < \. Observe that, for ah x G [0, A]'^, 
there exists j(x) G J"^ such that — "^/jl < ^fdi]/2. Therefore, 

\^^){x) - ^{x)\ < M^){iri) - 99(jr/)| + |$(v?)(j7?) - H^){x)\ 

+ W{iv)-'f{x)\ 

< r//2 + 2|jr? - 

<7?(\/d+l/2), 

so that the second part of the claim follows by substituting the value of r?. 

It remains to estimate the cardinality of := <l>(Lipf'). To see that, first 
observe that k[Qi) can take at most A/r/ values. On the other hand, it follows 
from (40) that, given A;(j), k{}-\-5i) can assume at most three different values, 
for all \ <l <d. This implies that 

< ^3(A/r,+l)'*-l ^ A 2\fd + 1 g((v/d+i/2)A/£+l)'' 

^ ~ rj e 6 



< 



A2\/d+2 ((^+l)A/e)d 



e 6 

the last inequality following since 1 < A/(2e). □ 

We can now estimate the error incurred when using an Euler approxima- 
tion of some future value of the empirical density process, centered on its 
current value. 

Lemma 6. For A; G Z+, n G N and a G [0, 1], 

¥{Wi{aMk + cTF(il/4), Mfe+LanJ ) > K^a^) < p, 

where a = 1 — a, K = Kp + 1, with Kp being the Lipschitz constant of F on 
V{X) in the variational distance, and 

Proof. First, observe that the following control of the increments holds: 
(42) \\Mk+i-Mk\\<2/n. 



28 G. COMO AND F. FAGNANI 

Define w := [crn\ and £ := KAa^. Also, for ip G Lipf^, define 

0<i<j 

for j = 0, . . . , w, and 

y(^) := (^Mk+^ _ (^1 _ Mfc - ^F(Mfc), - Z^f) . 
It follows from (42) that \\Mk+j - Mfc|| < 2j/n. Hence, 

|yM|=n"i Yl {F{Mk+j)-F{Mk),ip)- Y (MfcH_,-Mfe,v?) 

0<j<w 0<j<w 

<n"i ^ (||F(Mfc+,)-^(Mfc)ll + ll^fc+i-Mfcll)llvll 

0<j<w 



(43) 



^ n 

0<j<w 



<e/2, 

the last inequality following from the fact that \\ip\\ < A/2. Observe that, for 
all If G Lip;^(^), Zq'^^ = 0, while {Zj"^^ : < j < w} is a martingale. Moreover, 
(42) provides the following control on the increments: 

\Z^l\ - Zf\ < \{Mk+j+i - Mk+j,v) \ +n-'\{F{Mk+,) - Mk+j,^)\ 

(44) < \\Mk+j+i - Mk+j\\M\ +n~^\\F{Mk+j) - Mk+j\\M\ 
< 4n~^\\ip\\. 

Let 7i := 'Hs/12 ^ Lip^(^) be as in Lemma 5. By first applying the union 
bound, and then the Hoeffding-Azuma inequality ([2], Theorem 7.2.1) the 

probability of the event E := U/ig-h{I'^"''^^I — ^/4} can be estimated as fol- 
lows: 

(45) HE) < \n\n\Z\t^ > el 4) < m\ exp(-^^) . 
Now, Lemma 5 and (44) imply that 

Z(j~^^<3^\\^-h\\<3a±<'- 
for some h G l-Le/12- Hence, if E does not occur, then 

(46) |Z(r)| < min{|zW| + \Z^:^-'^\ ■.h^n}<'- 
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for every € Lipf^. By combining (43), (45) and (46), one gets 

nWi{Mk+n,.aMu + aF{Mk)) >e)= P(sup{z(r) + F^'^)} > e) 

<pfsup{z(r)}>^ 



<2|-H|exp 



and the claim follows upon substituting the expressions for w and e, and 
applying Lemma 5. □ 

We are now ready to prove point (b) of Theorem 1. Let L : = Lp — 1 and 
K = Kp + 1, where Lp and Kp are the global Lipschitz constants of F 
on V{X) in the Kantorovich-Wasserstein distance, and in the variational 
distance, respectively. Let us fix some e > 0, r > 0, and introduce the quan- 
tities 

Without any loss of generality, let us assume that a g]0,1], and put a = 
1 — 0". Further, let p be as in (41), and define 

(47) ao = e~^^V2, a^+i = (1 + (TL)ai + |KAo■^ iGZ+. 
Solving the iterative equation above, one obtains the estimate 

(48) «. = (! + aLY («o + ^ j - ^ < e'^'' [a, + 

For i £ Z"^, consider the random variable F" := Wi{Miw , fi^^i) , and the events 
Ai := {F" > ai}, Bi := Uo<j<i^i- shall prove by induction that 

(49) ¥{Bi) <{i + l)p 

for all i G Z4.. First, it follows from Lemma 3 that (49) holds with i = 0, 
for sufficiently small e and sufficiently large n. Then, for any nonnegative 
integer i, consider the intermediate measures 

A := aMwi + aF{Mwi), v := a/i^j + aF{[i„i). 

From the duality formula (39), and Lemma 2, one has 

(50) Wl(A,l/)<(^T + cJL^)F^ = (l + ^TL)^^ 

Furthermore, since {/ij} is a solution of the ODE (12), it follows from (15), 
and the estimate P^i(/U, < A/2\\fi — u\\, 

(51) \\pt-Pa^\\<2it-s), Wl{pt,Ps)<A{t-s) 



30 G. COMO AND F. FAGNANI 

for all t > s. From the duality formula (39), the fact that {fJ-t} solves the 
ODE (12), and (51), one gets the estimate 

W^i(i^,/i^(i+i)) = sup{(/i^(i+i),(/9) - (i/, : (/9 G Lipf } 

J ai 

(52) \\^lt-^lai\\dt 

ra{i+l) 
<AK / {t-ai)dt 

J ai 

= AKcjV2. 

From the triangle inequality, (50) and (52), one finds that 

rr+i < T^i(M^(,+i), A) + Wi{\, u) + Wl{u,^l,^^,+^)) 

(53) 

< H^i(M^(,+i), A) + K/\a^/2 + (1 + aL)V^. 

Therefore, (53), the inductive hypothesis (49), (47) and Lemma 6, imply 
that 

P(5,+i) = P(sfnAi+i) + P(Si) 

< ¥{Bf n {T^i(M^(,+i), A) > a,+i - KA,jV2 - (1 + aL)a^]) 
+ {i + ^)p 

< P(P^i(M^(,+i), A) > KAa2) + + l)p 

< ii + 2)p. 

Hence, (49) holds for all i G Z+. 

Observe that, if iw — w/2 < k <iw + w/2, then 

ni l^ai) 

(54) 

< + Aa. 

Now, recall the definition of <^(t) given in Lemma 4, and consider the events 
C := {?(r) < |nr} and D := {sup{|t - T^(t)\-te [0,r]} < e/(4A)}. Observe 
that C implies that, for all t <t, 



(55) tit) :-- 



<t) ^ 1 



[an\ 2 



3rn/2 1 2r 

< — + -<—. 

an — 1 2 a 
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It follows from (54), (51), (48) and (55), that, if the event -BfaraJ n n C 
occurs, then, for all t € [0,r], the following estimate holds: 

Wi{f,'l,f,t) = Wi{M,^t),f,t) 

< W^i(M,(i),/XT,,,,) + Tyi(/iT,(,),/Xt) 

<rrw + Aa + A|t-r,(i)| 

< a,^t) + Ad + e/4 
<e'^^^W + Acj + e/4 

<e2^"( ao + ^^ j +Aa + e/4 

where the last equality follows by substituting the expressions for a and ao- 
For sufficiently small e and large n, Lemma 4 implies that P(C n D) > 1 — p. 
Therefore, using (49), one gets that 

P(sup{M/i(/ir, f^t):te [0, r]} > e) < P(i?,(,)) + P(C^ U D') < (2r/a + 2)p, 

from which point (b) of Theorem 1 follows. 
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